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ABSTRACT 

We present an extension of the spectral deconvolution method ([Sparks fc Fordll2002l ) 
£f) ■ to achieve very high contrast at small inner working radii. We apply the method to the 

f^*) ' specific case of ground based adaptive optics fed integral field spectroscopy (without a 

C > . coronagraph). Utilising the wavelength dependence of the Airy and speckle patterns, 

we make an accurate estimate of the PSF that can be scaled and subtracted from the 
^ (— i data cube. The residual noise in the resulting spectra is very close to the photon noise 

O from the starlight halo. We utilise the technique to extract a very high SNR H & K 

I ' band spectrum of AB Dor C, the low mass companion to AB Dor A. By effectively 

eliminating all contamination from AB Dor A, the extracted spectrum retains both 
> ■ continuum and spectral features. The achieved la contrast is 9 mag at 0'.'2, 11 mag at 

^ | (y/5, in 20 mins exposure time, at an effective spectral bandwidth of 5.5 nm, proving 

that the method is applicable even in low Strchl regimes. 

The spectral deconvolution method clearly demonstrates the efficacy of image 
sheer based IFUs in achieving very high contrast imaging spectroscopy at small angular 
separations, validating their use as high contrast spectrographs/imagers for extreme 



C3 ' adaptive optics systems. 
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1 INTRODUCTION 

Direct imaging and spectroscopy of extra-solar planets is a 
research area that has attracted significant attention in re- 
cent years. Several specialised instruments are now being 
designed and built to achieve this goal within the next few 
years (e.g. VLT Plan et Finder (SPHERE) - iBeuzit et all 



(2006 ), Gemini GPI -ISoummer et ail l|2006l ). Keck XAOPI 
- iMacintosh et all i|2004l )). The major difficulty in detecting 
exo-planets is the extreme contrast between the star (typ- 
ically a main-sequence M-G dwarf) and the planet (typi- 
cally less than 1" away). Even in the most favourable con - 
ditions, theoretical computations (e.g. iBurrows et alj|2004l ) 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, under ESO programme 276.C-5013. 



have shown that the brightness ratio between the parent star 
and the exo-planet is expected to be at least 15 magnitudes 
at 1". Such contrasts are significantly beyond normal imag- 
ing techniques even with the best adaptive optics systems 
available today. If we are to attain this ambitious goal, novel 
ways to maximise the achieved contrast must be devised. 

Even in high-Strehl AO systems, long exposures do not 
provide a smooth halo around the central star, as one would 
expect from a time average of the speckles caused by at- 
mospheric turbulence. In fact, due to the presence of super- 
speckles, imaging contra sts are much worse than the theoret- 
ical photon noise limit (|Racine et al.lll999l). Super-sp eckles 
(also called quasi-static speckles Marois et alj ((2006)) are 
long lived speckles caused by imperfections in the light path, 
and have coherence times ranging from tens of seconds to 
several minutes. Any changes in the light path, for example 
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due to instrument flexure or changing telescope position, 
cause the super-speckles to vary, making it impossible to 
remove them with static calibrations obtained at a differ- 
ent time (e.g. PSF star observations will show a completely 
different super-speckle pattern). 

The contrast limitations posed by super-speckles can 
only be overcome through simultaneous observations. Sev- 
eral techniques have been proposed t o achieve this (e.g . 
iGuvonl 12004 ; iMarois et all 12004 120061 : iRen fc Wanel 12009 ) 
including the "spectral difference imager" (SDI) tech nique 
develop ed by one of the a uthors l|Lenzen et al.l [20041 ') and 
others l)Marois et al.ll2005h . The technique utilises the fact 
that the planetary spectrum has sharp, deep absorption fea- 
tures (e.g. CH4 in the H band) to dramatically improve 
the contrast with respect to scattered starlight, which has a 
smooth spectrum. By taking simultaneous images in and out 
of the methane absorption band, and scaling and subtract- 
ing the two images, one effectively removes the starlight, 
while leaving the planetary light intact. The subtraction re- 
moves residual starlight, whether in the form of rapidly time- 
varying speckles, slowly varying super-speckles, or any other 
form of scattered light within the telescope & instrument. In 
fact, the SDI technique uses an ingenious double-difference, 
so as to further reduce any residuals arising from syste matic 
errors due to differing optical paths. iBiller et all (l2006t ) have 
demonstrated the success of SDI in finding cool companions 
with the NACO instrument on the ESO-VLT. 

SDI is however limited in application as it relies on an 
intrinsic feature of the companion spectrum; the CH4 ab- 
sorption feature only found in objects with T c fr < 1200 K. 
In addition, there is a need for follow-up spectroscopy of the 
candidate object detected - both to confirm its nature (often 
via common proper motion with the parent), and to charac- 
teri se it in detail. In this paper, we describe an extension of 
the lSparks fc Fordl (|2002i ) spectral deconvolution (SD) tech- 
nique to achieve high-contrast imaging (and spectroscopy) 
with an integral field spectrograph (IFS). This technique 
does not rely on any feature of the target object's spectrum 
and is therefore applicable to any high contrast imaging ap- 
plication. In addition, the use of an IFS provides a complete 
dataset in one go, and with maximum signal-to-noise (no slit 
losses). Furthermore, complete 2D spatial information allows 
the continuum to be correctly measured, in contrast to ex- 
tremely narrow slits needed for A.O. assisted spectroscopy. 
No a-priori information about the spatial location of the 
companion is needed either, an inherent advantage of using 
an IFS. 

The next two section describes the ISparks fc Fordl 
i|2002l ) method, as well as proposed extensions, including 
using this technique at small inner working radii (within 
the bifurcation radius, see section [2~l2)) , and in section [3] we 
provide a demonstration of the achievable performance in 
this regime, even without a coronagraph. Section ^. 4| presents 
the results and section [4] elaborates on the future prospects 
and the discovery potential for high contrast imaging spec- 
troscopy with image sheer based integral field spectrographs, 
utilising this technique. 



2 CONCEPT 

Several authors (e g. ISparks fc Fordll2002l ; iFusco et al.ll2005l ; 
iBerton et alj|2006h have suggested that the wealth of spec- 
tral information available in an IFS data cube can be utilised 
to remove scattered starlight and identify the presence of 
a close-in companion, and extract its spectrum with en- 
hanced signal-to-noise (SNR), thus maximising contrast. 
Our scheme is an exte nsion of the spectral deconvolution 
technique proposed by ISparks fc Fordl (|2002T ), adapted to 
ground based observing with both extreme adaptive optics 
and general purpose adaptive optics systems, and further 
optimised to maximise SNR. Application of the concept to 
ground based IFS fed by general purpose A.O. systems is 
very challenging, as the Strehl ratio is typically only ~0.5, 
and varies strongly with wavelength. 

The central tenet of the SD technique is that the posi- 
tion of almost all features arising from the parent star (Airy 
pattern, speckles, etc.) scale smoothly and slowly with wave- 
length (e.g. the first Airy null is always 1.22 A/D from the 
star), whereas a physical object's location is independent of 
wavelength. Thus, stepping through the wavelength axis of 
an IFS data cube at a fixed spatial location (e.g. the location 
of a companion), one would see a modulation arising from 
maxima and minima of the Airy pattern p assing through 
the c ompanion's location (see figure 24 in ISparks fc Fordl 
2002). These modulations would drown out any expected 
signal from the companion. If however, we use the wave- 
length dependence to advantage, we can effectively subtract 
all wavelength dependent artifacts in the stellar point spread 
function (PSF), thus unmasking the presence of real physical 
objects. 

The technique is best illustrated via the removal of the 
Airy pattern from the IFS data cube, although it applies 
equally well to any achromatic aberration in the wavefront. 
Each wavelength slice of the IFS data cube is scaled radi- 
ally, so as to exactly compensate the wavelength dependence 
of the Airy pattern. Nulls and maxima in the Airy pat- 
tern, which used to be located at different radial distances 
in the different wavelength slices now line up perfectly at 
all wavelengths in t he scaled data cube (e.g. figure 27 in 
ISparks fc~F ord 2002). The companion's light now traces a 
dia gonal line through the data cube. 

ISparks fc Fordl ([2002) fit a low order polynomial to ev- 
ery spatial pixel (spaxel □) of the data cube, while rejecting 
outliers. The result is a very high SNR estimate of the Airy 
pattern (and all wavelength dependent speckles) at each 
wavelength of the data cube. The faint companion represents 
a high frequency modulation and is thus excluded from the 
low order polynomial fit (as it represents a tilted line through 
the scaled data cube, it is present in any given spaxel over 
only a small range of wavelengths) . The light from the par- 
ent star is then subtracted from each scaled slice of the data 
cube, and the result scaled back to the original grid, to form 
the result cube. 

The actual implementation of the concept for ground 
based A.O. assisted integral field spectroscopic data (with- 
out a coronagraph) has to account for low Strehl that varies 

1 The term spaxel describes all pixels spanning the entire wave- 
length range of an IFS data cube that correspond to a single 
spatial pixel 



Very high contrast IFU spectroscopy 3 

2.44 VD r V^i 



< *■ ► 




► e=l e>l 

'o 



Figure 1. The left hand panel is a schematic of a slice through the IFS data cube, in the x — A plane, with the x axis along the line 
joining the centre of the star with the centre of the faint companion. The dotted line represents the peak of the stellar light profile (as 
a function of lambda). The size of the companion, quantified as 2.44A/D, increases with wavelength. The middle panel shows the data 
cube scaled inversely with wavelength, centred on the star. The faint companion's separation is such that e = 1. The collapsed cube will 
be contaminated by companion light at all r values. The right hand panel illustrates the situation for epsilon > 1, where the hatched 
areas represent parts of the data cube not contaminated by light from the companion. It is thus possible to create a collapsed image free 
from companion light. 



substantially with wavelength. Furthermore, we have ex- 
tended the original concept in two ways: (a) to extract spec- 
tra from objects that lie close to the parent star and (b) op- 
timise the signal to noise ratio of the extracted spectra. We 
elaborate on these extensions in the following sub-sections. 

2.1 Application to faint companions at small 
separations 

The SD technique utilises the wavelength dependence of PSF 
features to obtain a high signal-to-noise estimate of the par- 
ent star PSF, while rejecting the light from the faint com- 
panion. We first derive the limiting distance for which the 
technique is applicable in its original proposed form, and 
then explain how we deal with close-in companions. 



2.2 Bifurcation point for characterisation analysis 

Consider a companion object located at radial distance r 
(expressed in angular units) from the primary object, im- 
aged with a telescope of diameter D. The distance to the 
first null of the Airy pattern is given by the usual formula- 
tion 

Go = 1.22^ (1) 

where Ao is the shortest wavelength in the IFS data cube. 
Defining the object extent as equal to 2€>o, and the wave- 
length range of the IFS as extending from Ao (shortest) to 
Ai (longest), we obtain an expression for the movement of 
the companion's centre in the scaled data cube as 

Ar = r — r — — r — — (2) 
Ai Ai 



Noting that the extent of the object stays constant in the 
scaled data cube, we can then express the bifurcation point 

as 



Ar = r^=2exl.22^ 
Ai D 



(3) 



where e is a factor slightly greater than one. The parameter 
e is explained in figure [T] 

For tsl, the scaled data cube is entirely contaminated 
by light from the companion, but only just so. For e > 1, 
a clean collapsed image devoid of flux from the companion 
object can be made by collapsing a small fraction of wave- 
length channels at either end of the data cube. This is only 
necessary in the near vicinity of the faint companion, as 
shown by the hatched areas in the right hand panel of figure 

m 

The amount by which e must exceed unity depends on 
the SNR of the observation, as sufficient number of wave- 
length channels must be used to form the high SNR PSF 
estimate. Table [1] lists values of the bifurcation point for 
typical near-infrared band passes for an 8 meter telescope, 
for a couple of typical values of e. It is obvious from the table 
that extended wavelength coverage by the IFS is crucial for 
removing the flux from a close-in companion in the scaled 
data cube. 



2.3 Extension to close-in faint companions 

We have extended the SD technique to close-in faint com- 
panions by using an iterative technique to remove the faint 
companion's light from the PSF estimate made using the 
scaled data cube. Residual light from the faint companion in 
the estimated PSF will reduce the flux level of the compan- 
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Table 1. Bifurcation point in the data reduction as a function 
of bandpass of observations and factor e. r is the value of the 
bifurcation radius in milli-arcseconds. 



Band 


Ai 


A 2 


c 


r 




/an 






mas 


H 


1.45 


1.8 


1.1 


516 


K 


1.95 


2.45 


1.1 


661 


H+K 


1.45 


2.45 


1.1 


246 


H 


1.45 


1.8 


1.2 


563 


K 


1.95 


2.45 


1.2 


721 


H+K 


1.45 


2.45 


1.2 


268 



ion in the final result cube. Furthermore, due to the wave- 
length scaling process, any residual light will affect spectral 
channels differently, resulting in an error in the continuum 
slope of the derived faint companion spectrum. Section 
shows that both the total flux and the continuum slope were 
correctly measured for the observations reported here, thus 
proving the efficacy of the proposed extension to the SD 
technique. 

The following assumes that the position of the close-in 
faint companion is already known (section l3.4l mentions how 
this might be derived from the IFS data cube itself). The 
companion is also assumed to be an unresolved point source, 
with a full width at half maximum (FWHM) roughly equal 
to the diffraction limited FWHM for the telescope aperture 
in question. 

Knowing both the position and the size of the source, 
we fit for the companion amplitude and background level 
in individual spectral channels of the observed data, prior 
to scaling. Subtracting an ideal Airy pattern centred on the 
primary improves the robustness of the fit, as does averag- 
ing ~10 spectral channels. Scatter in the fit parameters is 
further reduced by describing each fit parameter with a low 
order polynomial as a function of wavelength. The goal is 
to remove most of the flux from the faint companion, the 
quality of the fit need not be very good. The SD technique 
is then applied to the data cube from which the faint com- 
panion light has been removed, to first order. 

The result data cube will contain residual companion 
flux due to inaccuracies in the fit. However, as most of the 
PSF features will have been eliminated by the spectral de- 
convolution process, we can use the result data cube as a 
starting point for a second iteration, fitting the faint com- 
panion more accurately than in the initial step. After sub- 
tracting the improved estimate of the faint companion light 
from the observed data cube, we apply the SD method again 
(second iteration) to obtain an improved result cube. Fur- 
ther iterations are carried out until convergence is achieved, 
and the extracted spectrum of the faint companion does 
not change. Fit parameters such as companion position and 
width, fixed at first, can be free parameters in subsequent 
iterations to improve the quality of the fit. We applied this 
iterative technique to the data presented in the next section, 
and found that only two iterations were required to achieve 
convergence. 

The key to successfully applying the iterative procedure 
is to produce a reasonable, if noisy, estimate of the faint 
companion flux in the first iteration. This can be achieved 



by subtracting an ideal Airy pattern centred on the primary, 
or by subtracting data taken with the same instrument, but 
with a different rotator angle. In either case, the purpose 
is to estimate and remove the azimuthally symmetric com- 
ponent of the PSF, so as to enhance the contrast at the 
location of the faint companion. This can also be achieved 
by fitting a radial profile to the primary's light distribution, 
and subtracting it from the observed data cube. The radial 
profile can also be used to model the flux distribution from 
the faint companion, resulting in a better subtraction of the 
companion flux in the data processing. 

For faint companions located further than the bifurca- 
tion point, applying the iterative procedure outlined above 
can still improve the SNR of the final extracted spectra, as 
it reduces the contamination from the faint companion light 
in the PSF estimate. 

2.4 Optimising the SNR of the extracted spectra 

Observations carried out with a general purpose A.O. sys- 
tem suffer from the added complication of moderate Strehl 
that varies substantially with wavelength. Consequently, the 
FWHM of the PSF does not scale with A as expected for 
a purely diffraction limited system. Instead, with Strehl in- 
creasing with wavelength, the PSF is too broad at the short- 
est wavelengths, and too narrow at the longest wavelengths, 
a feature which persists in the scaled data cube. This re- 
sults in a radial distance and wavelength dependent flux 
residual after applying the SD technique. Consequently, the 
continuum of the faint companion is incorrectly determined. 
To rectify these errors caused by the varying Strehl, we fit 
and subtract a radial profile for the primary light from each 
spectral channel of the data cube prior to applying the SD 
technique. The data presented in section [3] have been so 
treated. 

A further improvement in SNR of the extracted spec- 
trum is obtained by scaling the PSF estimate cube back to 
the original observed pixel scale and subtracting it from the 
observed data cube, rather than scaling the observed data 
cube and performing the subtraction of the scaled cubes. 
Although the scaling is a linear operation, an interpolation 
is performed to place all scaled channels on a common grid, 
and this results in noise enhancement. As the observed data 
have much lower SNR (per spectral channel) than the PSF 
estimate formed by the SD technique, the noise enhance- 
ment is avoided almost completely by scaling and de-scaling 
the high SNR PSF estimate rather than the lower SNR ob- 
served data cube. 



3 DEMONSTRATION 

We obtained observations of the local young K-dwarf AB 
Doradus, which is known to have a M-dwarf com panion (AB 
Dor C) with AK=5 magnitudes only 0'.'2 away l|Close et al.l 
2005). AB Dor A has a K-band magnitude of 4.6, making it 
an ideal AO guide star. Details of the AB Dor C observations 
and data analysis are presented here as proof-of-concept of 
the proposed extensions to the SD technique, as described 
in previous sections. It is also the first time the SD method 
has been applied to a ground based A.O. assisted IFS data 
set, with moderate Strehl and no coronagraph present. As 
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Table 2. List of observations. All observations were made with 
the H+K grating. The 250 mas observations will not be discussed 
in this paper. 



Target 


Spaxel Size 
(mas) 


DIT 

00 


NDIT 


Total exposure 
00 


AB Dor a 


25 


5.0 


16 


1280 


AB Dor b 


25 


5.0 


16 


1280 


AB Dor c 


250 


0.83 


4 


6.64 


Teide l d 


250 


300 


1 


600 



a -16.5deg PA 
6 +16.5deg PA 

c These data were taken to bootstrap the spectral PSF from 250 
mas pixel scale to the 25 mas pixel scale, if required 
d Full designation is CI* Melotte 22 Teide 1, template of M8 
spectral type 



such, it also serves as an experimental verification of the 
SD method and its applicability to ground based data sets. 
We also present a high SNR H & K band spectrum of AB 
Dor C extracted using this technique. Further detailed anal- 
yses of the results of these observations are published in the 
companion paper bv lClose et all (|2007h (Paper II) . 



3.1 Observations 

Data were obtained on the nights of 24th a nd 25th Jan- 
uary 2006 with the S I NFONI instrument llThatte et all 
1 19981 ; lEisenhauer et al.l 120031 ; iBonnet et all |2004 ) at the 
Cassegrain focus of the ESO VLT-UT4 (Yepun). Tabled 
lists the observations, all with the H+K grating. DIT is the 
time per exposure, NDIT such exposures are averaged to- 
gether by the readout hardware, and the result is written to 
disk. The total exposure time is made up of many such aver- 
aged frames. In addition to the AB Dor system, observations 
were made of a star in the Pleiades cluster to serve as an 
M8 template of young age, and of several telluric standards, 
at both 25 and 250 mas plate scales. 

AB Dor A was used as th e guide star for th e SIN- 
FONI adaptive optics module i|Bonnet et al.l 120041 '). The 
mean Strehl ratio delivered was 0.37 on the 25th, and 0.35 
on the 26th, as reported by the SINFONI A.O. system itself. 
SINFONI was used in the 25 mas scale to ensure adequate 
sampling of the PSF. To enable removal of bad pixels, and 
to sample SINFONI's rectangular spaxels, a dither pattern 
with offsets of exactly ±0.5 spaxels was used. The H+K 
grating was used to maximise the wavelength coverage, re- 
sulting in a spectral resolving power of ~1500, and instan- 
taneous wavelength coverage of 1.4^tm- 2.45^im. No blank 
sky observations were taken for AB Dor at the 25 mas pixel 
scale, as it was clear that the data would be dominated by 
photon/speckle noise from AB Dor A. The J band was used 
for acquisition to avoid saturation of the detector array. See- 
ing during the observations (as reported by the Atmospheric 
Site Monitor) varied between 0'.'58 and 0'.'74, with a mean of 
/ .'64 on the 25th, and between 0'.'61 and / .'86 with a mean 
value of 0769 on the 26th. 



3.2 Data reduction 

Basic data reduction followed th e standard SINFONI d ata 
reduction procedure, outlined in Schrciber et all (|2004l ) & 
lAbuter et all {2006), with extensions to deal with the 2K 
camera upgrade of February 2005 jEisenhauer et al1l2002l ). 
The jitter pattern was corrected by combining data cubes, 
one for each on-source exposure recorded, creating a single 
mosaicked cube with a total observing time of 20 minutes 
for the data of the 25th, and a second one (with different 
rotator angle) for the 26th. Accurate centroiding of AB Dor 
A in the resulting data cube showed a residual variation 
with wavelength, implying that the standard differential at- 
mospheric refraction correction employed by the pipeline re- 
duction was not accurate enough. Consequently, we disabled 
the differential atmospheric refraction correction in the stan- 
dard reduction, and did a post-correction of the AB Dor A 
centroid movement along the wavelength axis of the data 
cube, using a cubic polynomial function. The resulting im- 
proved model for differential atmospheric refraction will be 
presented elsewhere. 

It should be noted that even after disabling the dif- 
ferential atmospheric refraction correction, there are still 
three interpolations inherent in the production of a SIN- 
FONI data cube, two of which use a nearest neighbour in- 
terpolation algorithm with a tanh kernel, and one uses a 
polynomial interpolation algorithm. One interpolation is re- 
quired to "straighten" the spectra, a second does re-gridding 
onto an uniform grid in wavelength, and a third corrects 
for sub-pixel shifts between slices. We have not evaluated 
the impact of these interpolations in a quantitative manner, 
but it is likely that they do impact the maximum contrast 
achieved by these observations. 

Implementation of the SD technique (and its exten- 
sions) used both IDL and GIPSY data analysis packages, in 
addition to the routines fro m the SINFONI pipeline. GIPSY 
l|Vogelaar fc Terlouwll200"lh was used primarily for visualisa- 
tion of the data cubes with sliceview, while most compu- 
tationally intensive tasks were performed in IDL. 

3.3 The SD technique and extensions 

AB Dor C is located 0'.'2 from AB Dor A, at P.A. 181° at 
the epoch of t hese observations, ba sed on orbit parameters 
determined by iNielsen et al.l l|2005h . Using equation [3] we 
note that AB Dor C is closer than the bifurcation radius, 
thus requiring an iterative procedure to identify the com- 
panion and extract its spectrum. 

The technique described in previous sections was ap- 
plied to the AB Dor data cubes to extract a spectrum of the 
companion AB Dor C. As AB Dor C is clearly visible in the 
raw data cube (see figure [2]), there was no problem in de- 
termining the companion location. However, AB Dor C did 
overlap with a bright spot in the diffraction pattern in the 
data taken on the 25th. The pattern of four bright spots is 
caused by the telescope secondary mount diffraction pattern 
superposed on the Airy pattern. Consequently, accurate de- 
termination of the companion location did require a second 
iteration of the SD technique, as explained in section 12.31 
above. Due to the strong dip in atmospheric transmission 
at the edges of both H and K bands, the processing was 
confined to wavelength ranges 1.457 to 1.80 and 1.95 to 
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5 h 28 m 44.95 s 44.90' 



R.A. (J2000) 

Figure 2. Illustration showing the efficacy of the SD technique at removing both PSF artifacts and super-speckles from the SINFONI 
IFS data cube for the AB Dor system. The left frame shows one wavelength slice of the observed data cube at 2.2 /im. Note that the 
entire vertical extent of the image is only 0"9. The colour table is logarithmic (min 10 , max 10 4 ). The contours are logarithmic, from 
0.9 to 2.3 in steps of 0.1, and from 2.3 to 4.0 in steps of 0.2. The middle frame shows the same data, but with a radial profile fitted 
and subtracted, so as to highlight the PSF imperfections. The super-speckles are easily confused with real sources in this narrow-band 
image. The colour table is now linear (min -10, max 25), with contours from -12.5 to 32.5 in steps of 5, and from 32.5 to 150 in steps 
of 20. The four-fold symmetry of the Airy pattern arises from the superposition of the diffraction spikes from the secondary support 
structure on the Airy rings. The right hand frame shows the same wavelength slice of the data cube, after applying the SD technique 
iteratively. Colour table and contours are the same as for the middle frame. Super-speckles are completely absent at the lowest contour 
level of ±2.5, corresponding to a lcr error of <1 unit. 



2.45 /xm respectively for the H and K bands. A total of two 
iterations were needed for the SD technique to yield a high 
SNR PSF estimate free from companion flux. 

3.4 Results 

Figure [2] shows one wavelength slice of the SINFONI IFS 
data cube for the AB Dor system before and after the ap- 
plication of the SD technique (and its extensions) described 
above. The separation of AB Dor C from the parent was 
only C/.'2 at this epoch, while the flux contrast was Am = 
4.90. The left hand frame shows the raw data, while the cen- 
tral frame has a smooth radial profile fitted and subtracted 
from the raw data to highlight the PSF artifacts, and in 
particular, the super-speckles. The right hand frame shows 
the result of applying the technique described in section [2] 
We can successfully remove both Airy pattern residuals, as 
well as super-speckles. Note the four-fold symmetry of the 
Airy pattern, a result of the superposition of the Airy ring 
with the diffraction from the spider arms of the secondary. 
Unfortunately, AB Dor C is located at the same orientation 
as one of the Airy peaks. We prefer to use this data set as 
it has significantly better Strehl, although the other data 
set (rotated by 33°) places AB Dor C at a more favourable 
location relative to the Airy peaks. Stepping through the 
data cube along the wavelength axis, all PSF imperfections 
scale linearly with wavelength, Q while the real object stays 

2 mpeg movies of the raw, radial profile sub- 
tracted and SD result data cubes can be viewed at 
http : //www-astro .physics . ox . ac.uk/ ~thatte/abdorc 



at a fixed location. The data do demonstrate that the SD 
technique can effectively separate companion and PSF con- 
tributions, even if they are spatially coincident in a single 
wavelength channel. The removal of the super-speckles is of 
particular value, as their shape mimics point sources in the 
field in any narrowband image. 

3.4-1 Extracted spectra of AB Dor C 

Figures [3] and [4] shows the normalised, extracted spectra of 
AB Dor C, using the SD technique. The data have been 
smoothed with a boxcar of width 3 pixels, corresponding to 
1.5 nm. The SNR of the K band spectrum exceeds 40, as de- 
rived from the contrast values discussed in !3.4.2l The quality 
of the spectrum is limited by systematic errors in the division 
of telluric features, as no sky observations were performed. 
The N ACO spectrum previously obtained by IClose et al.l 
l|2005h at a tighter / .'l5 separation is also shown for com- 
parison. Note that the SD technique is able to reproduce not 
only the spectral features but the continuum slope as well, 
in contrast to A.O. long-slit measurements with varying slit 
filling factors. The output data cube, after applying the SD 
technique, is used to measure the position, flux, and spec- 
tral type of AB Dor C. At the epoch of these observations, 
AB Dor C was located 201.7±10 milli-arcseconds from AB 
Dor A, at P.A. 180.78°. The K magnitude was 4.90 mag 
fainter than AB Dor A, which, combined with the 2MASS 
measurement for AB Dor A yields K=9.59 for AB Dor C. A 
spectral type of M5.5±0.5 can be derived by comparing the 
observed spectrum with young and old templates. A detailed 
analysis of the properties of AB Dor C (including its spec- 
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Figure 3. The K band extracted spectrum of AB Dor C, with 
some of the prominent stellar features marked. The K band spec- 
trum obtai ned by NACO A.O. long slit spectroscopy, at 0'.'150 
separation ^Nielsen et alj l2005f) is also shown for comparison. 
Note that our spectrum also correctly recovers the continuum 
slope, vital for deriving the spectral type of this young object 
(see Paper II for details of the spectral classification). The "emis- 
sion" feature just shortward of the first CO bandhead is a residual 
telluric feature that was not correctly subtracted out. 



tral type determination) is presented in Paper II. Note that 
both the companion flux and its spectral continuum slope 
are correctly derived using the iterative method described 
in section [2731 above. 



3.4-2 Achieved Contrast 

The remarkable success of the SD technique at removing 
light contamination resulting from PSF imperfections is il- 
lustrated in figure [5] The plot shows the resulting radial pro- 
file for a 20 min exposure of the AB Dor system observed 
with SINFONI, with a spectral width of 5.5 nm, using the 
SD technique. AB Dor C is clearly visible as a bump in the 
radial profile at 0'.'2. The curve shows the contrast achieved 
(la) is ~ 9 mag at 0'.'2, and 11 mag at 0'.'5, without us- 
ing a coronagraph. Note that the effective spectral resolu- 
tion (R=400) is higher than that achi eved by SDI (11=50) 
jLenzen et all 12004 iBiller et al.l |2006| ). The achieved con- 
trast is still ~1 mag from the photon noise limit, although 
we suspect that some residual noise is attributable to non- 
optimal interpolations in the data reduction procedure, and 
can be improved upon. The ultimate limitations of the tech- 
nique will be addressed in a separate paper. 



3.5 Faint object detection 

AB Dor C is only 5 mag fainter than AB Dor A, and could 
thus be easily detected in the raw data cube. However, we 
note that the SD technique can also be used for detection 
of faint companions. Any feature whose radial distance from 
the primary does not scale with wavelength generates a ra- 
dial streak in a collapsed image created by averaging the 
scaled data cube along the spectral dimension. Subtract- 
ing the collapsed image (scaled back appropriately for each 



Figure 4. The H band extracted spectrum of AB Dor C, with 
some of the stellar features marked. Also marked are the location 
of the Brackett series hydrogen recombination lines, as these were 
present in absorption in the late B telluric standard. Although ev- 
ery effort was taken to accurately fit these lines, some residuals 
remain. As no sky observations were made, some of the bright- 
est night sky OH emission lines are also present in the derived 
spectrum. The H and K spectra were obtained simultaneously, 
by applying the SD technique (and its extensions) to the data 
cube covering both H & K bands. They are shown separately for 
clarity. 



VLT/SINFONI: PSF Scaling and Differencin 





0.0 0.2 0.4 0.6 0.8 

Distance from star (arcsecs) 

Figure 5. The plot shows three curves - the top curve is the 
radial profile of AB Dor A. The bottom curve is the square root 
of the top curve, so it represents the photon noise limit. The 
middle curve is the standard deviation in the SD result frame, so 
it is a measure of the residual noise. The horizontal lines at 6, 8 
and 10 magnitudes are for reference. The large bump at 0'.'2 in 
the SD result frame is due to AB Dor C. The achieved contrast 
is 9 mag at 0"2, and 11 mag at 0'.'5, in 20 mins exp. time at 
R e ff=400. As no coronagraph is used, very high contrasts can be 
obtained at small inner working radii. 
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wavelength) from the original data cube results in a dimin- 
ished companion flux plus an inward negative streak at the 
shortest wavelengths, a diminished flux companion plus a 
symmetric negative radial streak at central wavelengths, and 
a weaker companion plus an outward negative streak at the 
longest wavelengths. The final collapsed image shows the 
characteristic pattern of a positive object sitting atop a long 
negative radial st reak centred on the ob ject ((see both "plan- 
ets" in fig. 26 of ISparks fc Fordll2002l )). This characteristic 
pattern can be used to detect faint companions. Further- 
more, if the primary's radial light profile is fitted and sub- 
tracted from each wavelength slice of the data cube prior to 
applying the SD technique, the contrast of the faint com- 
panion can be further enhanced. 

If the technique is applied to data cubes taken at two 
different rotator angles and subtracted from each other prior 
to step 1, only features that are stationary on the sky will 
remain. In addition to showing the characteristic pattern 
described above, they will also show a positive and negative 
image of true celestial objects, separated by the difference 
in rotator angles. The resulting pattern can be easily used 
to identify potential faint companions. 



4 CONCLUSIONS 

In this paper we h ave introduced extensi ons to the SD tech- 
nique proposed bv lSparks fc Fordl (|2002h for achieving high 
contrast imaging spectroscopy with an AO-fed integral field 
spectrograph at small inner working radii. Applying this 
technique to real data, we have shown it to provide very high 
contrast spectra of faint companions very close to bright 
stars. The achieved contrast is substantially better than pre- 
viously achieved by other techniques, especially as no coron- 
agraph was used. The absence of a coronagraph removes lim- 
its on the smallest inner working radius, especially valuable 
as most exo-planet candidates are expected to lie at very 
small separations from the parent star. The SD technique 
holds great promise for direct imaging of exo-planets, as it 
simultaneously detects and spectrally characterises any faint 
companion, thus removing the need for expensive and time- 
consuming follow-up observations, either to detect common 
proper motion or to obtain a spectrum of the faint compan- 
ion. In addition, the technique does not require any assump- 
tions about the companion's spectral characteristics (e.g. 
presence of a CH4 feature), and can therefore be applied to 
any high contrast application. 

Our demonstration of the efficacy of the SD tech- 
nique to obtain spectra of the close-in faint companion AB 
Dor C shows that image slicer based integral field spectro- 
graphs are capable of achieving very high co ntrasts, con- 
trary to the expectations o f several groups |Claudi et al.l 
12004 IVerinaud etakl 120061 ; | Berton et al. I I2006F who were 
concerned that large non-common-path errors would limit 
the contrast achievable with these s ystems. The SINFONI 
IFS design (|Tecza et al.l l200l I2003T I exclusively uses clas- 
sically polished flat glass mirrors in its image slicer, thus 
achieving very small non common path errors, as demon- 
strated here. Indeed, we have shown that an integral field 
spectrograph alone can provide a large fraction of the to- 
tal contrast requirement of an exo-planet direct detection 
instrument. 



Successful application of the SD technique requires large 
instantaneous wavelength coverage, strongly favouring im- 
age slicer based IFS designs. These naturally provide large 
simultaneous wavelength coverage, while it is rather difficult 
to achieve a large bandwidth in lenslet array based designs. 
The requirement to keep non common path errors to an ab- 
solute minimum strongly favours image slicer designs with 
classically polished glass slicing optics (not just slicing mir- 
rors, but the entire slic er optics), as implemen ted in the 
SINFONI spectrograph (|Tecza et al.ll200dl2003h . 
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